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Abstract. The current picture of the collapse and explosion of massive stars and the 
formation of neutron stars is reviewed. According to the favored scenario, however by 
no means proven and undisputed, neutrinos deposit the energy of the explosion in 
the stellar medium which surrounds the nascent neutron star. Observations, in par- 
ticular of Supernova 1987A, suggest that mixing processes play an important role in 
the expanding star, and multi-dimensional simulations show that these are linked to 
convective instabilities in the immediate vicinity of the neutron star. Convectively en- 
hanced energy transport inside the neutron star can have important consequences for 
the neutrino emission and thus the neutrino-heating mechanism. This also holds for a 
suppression of the neutrino interactions at nuclear densities. Multi-dimensional hydro- 
dynamics, general relativity, and a better understanding of the neutrino interactions in 
neutron star matter may be crucial to resolve the problem that state-of-the-art spher- 
ical models do not yield explosions even with a very accurate treatment of neutrino 
transport by solving the Boltzmann equation. 



1 Introduction 

Baade and Zwicky H were the first who speculated about a connection between 
supernova explosions and the origin of neutron stars. They recognized that stellar 
cores must become unstable because neutrons, being produced by captures of 
degenerate electrons on protons, define an energetically advantageous state. The 
gravitational binding energy liberated by the collapse of a stellar core could 
power a supernova explosion. More than thirty years later, Colgate and White 
explored this idea by performing numerical simulations pl| | . Their models showed 
that the prompt hydrodynamical shock, which forms at the moment of core 
bounce, is not able to reach the outer layers of the star due to severe energy losses 
by photointegration of iron nuclei. Therefore they suggested that the neutrinos, 
which arc emitted in huge numbers from the hot, collapsed core and carry away 
the binding energy of the forming neutron star, deposit a fraction of their energy 
in the stellar medium external to the shock. Later, more accurate simulations 
with an improved treatment of the input physics such as nuclear equation of state 
and neutrino transport |P, pTj|Jll] . [l^ . [l9p |j60t| confirmed the failure of the prompt 
shock, but could not find the neutrino-driven explosion imagined by Colgate and 
White. 

The discovery of weak neutral currents, whose existence has been predicted 
by the standard model of electro- weak interactions, had a major influence also on 
supernova models. It became clear that the neutrinos are strongly coupled to the 
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stellar plasma and can escape only on timescales of seconds. This implied that 
the neutrino luminosities did not become large enough that the rate of energy 
or momentum transfer by neutrino interactions with the nuclei in the stellar 
plasma ahead of the supernova shock could become dynamically important. In 
1982, however, Jim Wilson discovered (63) that the stalled prompt shock can be 
revived by neutrino heating in the dissociated postshock medium on a longer 
timescale. While immediately after shock breakout neutrino emission extracts 
energy from the shocked gas, energy transfer from neutrinos to the postshock 
medium is favored hundreds of milliseconds later. If the shock has expanded to a 
larger radius and the postshock temperature has decreased, energetic neutrinos, 
which stream up from deeper regions, can deposit a small fraction of their energy 
in a gain layer behind the shock @j. Although the efficiency is much lower than 
imagined by Colgate and White — typically only a few per cent of the neutrino 
energy are left in the stellar medium — the neutrino heating mechanism can 
yield the energy of about 10 51 erg for the explosion of a massive star like, e.g., 
Supernova 1987A. 

Although Wilson's simulations showed the principle viability of this mecha- 
nism and theoretical investigations enlightened the underlying physics jl|,[3]]£],|7j , 
later simulations revealed a strong sensitivity to the details of the post-bounce 
evolution of the collapsed stellar core |^,||,^3[^] . General relativity, the nuclear 
equation of state and corresponding properties of the nascent neutron star, the 
treatment of the neutrino transport and neutrino-matter interactions, and the 
structure of the collapsing star have important influence. In particular, later 
models of Mayle and Wilson |34],[35) produced sufficiently powerful explosions 
only in case of enhanced neutrino emission from the nascent neutron star due to 
neutron-finger instabilities. Whether such mixing processes take place, however, 
is not clear to date [jl5| . Alternatively, the hot neutron star might develop Ledoux 
convection J2^,[^] due to negative gradients of lepton number and/or entropy, 
a possibility which is suggested by the structure obtained in state-of-the-art 
spherical models for the neutrino-cooling phase of proto-neutron stars |50| , ^8| 
and by two-dimensional hydrodynamical simulations |p8[ , but is not generally 
accepted ]44| . 

Supernova 1987A has provided us with evidence for large-scale mixing pro- 
cesses, which involve the layers in the near vicinity of the neutron star where ra- 
dioactive nuclei are produced. Indeed, multi-dimensional simulations have shown 
the existence of hydrodynamical instabilities in the neutrino-heated region al- 
ready during the first second of the explosion 1 28 47],|32|,[53 20 . This convective 



overturn behind the shock was recognized to allow for an explosion even when 
spherical models do not explode. Also, it enhances the explosion energy and sets 
the timescale until the shock gains momentum. But still it is not clear whether 
it is sufficient and crucial for the neutrino-driven mechanism to be at work, or 
whether it is just an unavoidable side-effect jl5). The existence of these mixing 
processes, however, definitely means that the events leading to an explosion, 
the energetics of the explosion, the nucleosynthesis of radioactive elements dur- 
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ing the first second, and the composition and distribution of the ejecta can be 
understood quantitatively only by multi-dimensional simulations. 

Unfortunately, the 24 neutrinos from Supernova 1987A, which were recorded 
in the underground experiments of Kamioka, 1MB and Baksan Q , did not pro- 
vide enough statistical information to draw conclusions on the events which 
instigated the explosion of the star. Although these neutrinos confirm the basic 
picture of stellar core collapse and neutron star formation, they are not suitable 
to support the neutrino-driven mechanism. Therefore the latter may be consid- 
ered as the currently most favorable explanation for the explosion, but empirical 
evidence cannot be put forward and numerical models do not draw a clear and 
unambiguous picture. 

It cannot be excluded that the energy for supernova explosions of massive 
stars is provided by some other mechanism than neutrino heating, for example 
by magneto-hydrodynamical processes. However, we know that neutrinos with 
the expected characteristics are emitted from collapsed stellar cores, we know 
that these neutrinos carry away the gravitational binding energy of the nascent 
neutron star, we know that neutrino heating must occur behind the stalled shock 
some time after core bounce, we know that analytic studies and numerical sim- 
ulations find explosions for a suitable combination of conditions, we know that 
the energy transfer by neutrinos can be strong enough to account for "normal" 
explosions with a canonical energy of ~ 10 51 erg, and we know that the timescale 
of shock rejuvenation by neutrino heating determines the mass cut such that ac- 
cretion of the collapsed stellar core and later fallback lead to neutron star masses 
and supernova nucleosynthesis in rough agreement with observations [p3|j2^j67|| . 
Of course, the discrepant results of numerical simulations are unsatisfactory, 
and major problems are nagging: Why do the supposedly best and most ad- 
vanced spherical models not produce explosions? Can one trust current multi- 
dimensional simulations with their greatly simplified and approximate treatment 
of neutrino transport? Is rotation in neutrino-driven explosions sufficient to ex- 
plain the observed large asphericities and anisotropics in many supernovae |24j ? 
What is the reason for the kicks by which pulsars are accelerated to average 
velocities of several hundred km/s presumably during the supernova explosion? 
What powers hyperenergetic supernovae, which seem to release up to 50 times 
more kinetic energy than ordinary explosions of massive stars 1 30 ? 



This paper gives an overview over the major lines of research on the stan- 
dard explosion scenario of massive stars, and on neutron star formation, where 
progress has been achieved during the past few years or will be coming up in the 
near future. In Section 2 the physics of neutrino-driven explosions will be dis- 
cussed in some detail, and an analytic toy model will be described which allows 
one a deeper understanding of the requirements of neutrino-driven explosions. 
In Section 3 the current status of supernova modeling in spherical symmetry will 
be outlined, with particular focus on results with neutrino transport by solving 
the Boltzmann equation. In Section 4 the relevance of convective overturn in 
the neutrino-heated layer will be discussed. New results of two-dimensional sim- 
ulations will be presented, which take into account the nucleosynthesis during 
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Fig. 1. Sketch of the post-collapse stellar core during the neutrino heating and shock 
revival phase. R v is the neutrinosphere radius, R na the protoneutron star radius, f? g 
the gain radius outside of which net neutrino heating exceeds neutrino cooling, and 
R s is the shock radius. The shock expansion is impeded by mass infall at a rate M, 
but supported by convective energy transport from the region of strongest neutrino 
heating into the post-shock layer. Convection inside the nascent neutron star raises the 
neutrino luminosities. 

the explosion and follow the shock from the moment of its formation until it 
breaks out of the surface of the star. In Section 5, two-dimensional hydrody- 
namical models of the neutrino-cooling phase of newly formed neutron stars will 
be described. The effects of rotation and of a suppression of neutrino-nucleon 
interactions by nucleon correlations in the nuclear medium will be addressed. 
Section 6 will contain a summary and conclusions. 

2 The Explosion Mechanism 

The physics of neutrino-driven explosions is discussed, first on the level of basic 
considerations, then with the help of an analytic toy model, which allows one to 
study the competing effects that determine the destiny of the stalled supernova 
shock. 

2.1 Neutrino Heating 

Figure [l] displays a sketch of the neutrino cooling and heating regions outside the 
proto-neutron star at the center. The main processes of neutrino energy deposi- 
tion are the charged-current reactions v e + n — > p + e~ and v e + p — > n + e + || . 
With the neutrino luminosity, L, the average squared neutrino energy, (e 2 ), and 
the mean value of the cosine of the angle Q v of the direction of neutrino propa- 
gation relative to the radial direction, (fi) = (cos#j,), being defined as moments 
of the neutrino phase space distribution function f(r, t, /i, e) by integration over 
energies e and angles /i according to 

/+1 poo f p + 1 poo ~i ~! 

^ d M y o dee5 f-\J 1 d ^ J o dee3 f\ > ( 2 ) 
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where Y„ and Y p are the number fractions of free neutrons and protons (number 
densities divided by baryon number density), respectively. In the second equation 
Y n +Yp w 1, and equal luminosities and spectra for v e and v e were assumed. L v ^2 
denotes the total luminosity of v e plus v e in 10 52 erg/s, r 7 the radial position 
in 10 7 cm, (e^ 15 ) is measured in units of 15 MeV, and / = (/x)„ is very small 
in the opaque regime where the neutrinos are isotropic, adopts a value of about 
0.25 around the neutrinosphcre, and approaches unity for radially streaming 
neutrinos very far out. Note that the "flux factors" {\x) v determines the neutrino 
energy density at a radius r according to e u (r) — L„/(47rr 2 c(/i)„). 

Using this energy deposition rate, neglecting energy losses due to the re- 
emission of neutrinos, and assuming that the gravitational binding energy of 
a nucleon in the neutron star potential is (roughly) balanced by the sum of 
internal and nuclear recombination energies after accretion of the infalling matter 
through the shock, one can estimate (very approximately) the explosion energy 
to be of the order 
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Here AM is the heated mass, At the typical heating timescale, Ego the (net) 
total gravitational binding energy of the overlying, outward accelerated stellar 
layers, and E nnc the additional energy from explosive nucleosynthesis, which is a 
significant contribution of a few 10 50 erg only for progenitors with main sequence 
masses above 20 M Q , and which roughly compensates -EgbQ 

It is not easy to infer from Eq. the dependence of the explosion energy 
on the neutrino luminosity. On the one hand, 



£ox P cx QtVAt 



AM At cx L v At At, 



(6) 



The latter statement is supported by the following argument (S. Woosley, personal 
communication): Material with a specific gravitational binding energy # gra v which 
is equal to or larger than the nuclear energy release per gram in Si burning, e nuc ~ 
10 18 erg/g, is located interior to the radius where the temperatures can become high 
enough (T J> 5 x 10 9 K) for explosive nucleosynthesis of 56 Ni. ^From <£ gra v = GM/r J> 
10 18 erg/g one estimates a radius of r ^ 2 x 10 s cm, and from ^nr 3 aT 4 ~ 10 51 erg 
with I ^ 5 x 10 9 K one finds r <, 4 x 10 s cm. This means that the energy release 
from explosive nucleosynthesis is easily able to account for the gravitational binding 
energy of the burning material. 
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where V is the heated volume between gain radius and shock, and At the optical 
depth of the heating layer. On the other hand, 

At cx (el)R gPg cx (el)R g T* cx ^(e 2 ,) 3 / 2 . (7) 

Here p cx T 3 was assumed for the relation between density and temperature 
in the heating layer |j, and Q+(i? g ) = Q~(R g ) was used at the gain radius, 
where neutrino heating is balanced by neutrino cooling. The energy loss rate 
Q~ by neutrinos produced in capture reactions of nondegenerate electrons and 
positrons on nucleons, scales with T 6 . Combining Eqs. (|^) and (Q) yields 

£ cxp oc L^{el)^At oc L^At, (8) 
when (e 2 ) cx T 2 cx L]J 2 is used for black-body like emission. 

1/2 

If the expansion velocity were simply proportional to -E C xp, in which case 

— 1/2 

the time At for the shock to reach a given radius would be At cx -Ecxp , then 
i^cxp cx Lu @- However, when shock expansion sets in, most of the energy is 

1/2 

internal energy, but not kinetic energy, making the relation At cx i? C xp very 
questionable. The actual variation of At with the inverse of the neutrino lumi- 
nosity can be steeper. 



2.2 Requirements for Neutrino-Driven Explosions 



In order to get explosions by the delayed neutrino-heating mechanism, certain 
conditions need to be fulfilled. Expansion of the postshock region requires suf- 
ficiently large pressure gradients near the radius i? C ut of the developing mass 
cut. If one neglects self-gravity of the gas in this region and assumes the den- 
sity profile to be a power law, p(r) cx r~ n (which is well justified according to 
numerical simulations which yield a power law index of n ~ 3 §), one gets 



P(r) 



for the pressure in an atmosphere near hydrostatic equilibrium. 



Outward acceleration is therefore maintained as long as the following condition 
for the "critical" internal energy density e holds: 
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where use was made of the relation P = (7 — l)s. The numerical value was ob- 
tained for 7 = 4/3 and n = 3. This condition can be converted into a criterion for 
the entropy per baryon, s. Using the thermodynamical relation for the entropy 
density normalized to the baryon density rib, s = (e + P)/(ribT) — J2i ViYi where 



r)i (i 



p, e , e + ) are the particle chemical potentials divided by the temper- 



ature, and assuming completely disintegrated nuclei behind the shock so that 
the number fractions of free protons and neutrons are Y p = Y e and Y n — 1 — Y e , 
respectively, one gets 
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In this approximate expression a term with a factor Y e was dropped (its absolute 
value being usually less than 0.5 in the considered region), nucleons were assumed 
to obey Boltzmann statistics, and T is measured in MeV, Mi in units of M©, 
pg in 10 9 g/cm 3 , and r-j in 10 7 cm. Inserting typical numbers (T w 1.5 MeV, 
Y n 0.3, R cu t ~ 1.5 • 10 7 cm), one finds s > Ibk^/N, which gives an estimate 
of the entropy in the heating region when expansion is going to take place. 

Since the entropy and energy density in the postshock later are raised by 
neutrino energy deposition, the conditions of Eqs. (^) and (|Io|) imply require- 
ments on the neutrino emission of the proto-neutron star. These can be derived 
by the following considerations. A stalled shock is converted into a moving one 
only, when the neutrino heating is strong enough to increase the pressure behind 
the shock by a sufficient amount. From the Rankine-Hugoniot relations at the 
shock, a criterion can be deduced for the heating rate per unit mass, q v , behind 
the shock, which leads to a positive postshock velocity (u\ > 0) 

q > 2 Izl (ii) 

q " > /33(/3-l)( 7 -l) V R S ■ {U > 

Here (3 is the ratio of postshock to preshock density, (3 = pi/po, 7 the adiabatic 
index of the gas (assumed to be the same in front and behind the shock), and 
r\ defines the fraction of the shock radius R s where net heating by neutrino 
processes occurs: 77 = (R s — R g )/R s . uo is the preshock velocity, which is a 
fraction a (analytical and numerical calculations show that typically a ~ 1/ v2) 
of the free fall velocity: uq — ay / 2GM/r. Assuming a strong shock, one has 
= (7 + l)/(7 — 1), which becomes f3 = 7 for 7 = 4/3. With typical numbers, 
R s — 100km, M = Mi = lAf Q , and 77 w 0.4, one derives for the critical 
luminosity of v e plus v e : 

M 3/2 

L»M<i 5 ) > 4 -4 ^7T- (12) 

Since this discussion was very approximate, e.g., the reemission of neutrinos was 
ignored and properties depending on the structure of the collapsed stellar core 
were absorbed into free parameters, the analysis cannot yield a quantitatively 
meaningful value for the threshold luminosity. However, the existence of a lower 
bound on the neutrino luminosity as found in numerical simulations |^,^3|, is 
confirmed. Above this threshold value, neutrino heating of the gas behind the 
shock is strong enough to drive an expansion of the gain layer. 



2.3 Analytic Toy Model 

The discussion in the previous section was overly simplified. The behavior of the 
stagnant shock does not only depend on the neutrino heating in the gain region, 
but is also influenced by the energy loss and the settling of the cooling layer, 
and by the mass infall from the collapsing star star ahead of the shock. Only 
hydrodynamical simulations can determine the shock evolution in response to 
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Fig. 2. Left: Phase diagram for successful explosion or continued stellar collapse in 
dependence upon the rate at which gas falls into the shock and upon the v e plus v e 
luminosity of the nascent neutron star. The mass of the neutron star was assumed to 
be 1.25 Mq. Right: Shock positions as functions of time for different values of the v e 
plus i? e luminosity of the neutron star (according to the labels) and fixed value of the 
rate at which gas falls into the shock (marked by the dashed line in the left plot). For 
neutrino luminosities above the critical value (where the dashed line crosses the solid 
line in the left plot) explosions can occur. 



these different, partly competing effects. Analytic discussions, however, can help 
one understanding the significance of different effects and thus can supplement 
more detailed, but less transparent supercomputer calculations. 

The conditions in the region of neutrino heating and the details of the heating 
process were analysed by Bethe and Wilson || and Bethe j|,f|jg,[7| ■ Burrows and 
Goshy jl8) considered the post-bounce accretion phase of the supernova shock 
as a quasi steady-state situation, and thus replaced the time-dependent partial 
differential equations of hydrodynamics by a set of ordinary differential equations 
to determine the radial position of the standing shock by an eigenvalue analysis. 
While this approach captures an interesting aspect of the problem, it has serious 
weaknesses. The accretion flow between the shock and the neutron star does not 
need to be stationary, but neutrino-heated matter may stay in the gain region, 
or gas will pile up on the forming neutron star when neutrinos are unable to 
remove energy quickly enough for the gas to settle. In particular, when the 
shock accelerates outward, the transition from accretion to outflow cannot be 
described as a stationary situation. 

Since the gas falling into the stalled shock is strongly decelerated and the 
postshock velocities are smaller than the local sound speed, the structure of the 
collapsed stellar core is rather simple and can be well described by hydrostatic 
equilibrium. This allows for an approximate treatment, which is complementary 
to the approach taken by Burrows and Goshy []l8| : Integrating the stellar struc- 
ture over radius leads to conservation laws for the total mass and energy in the 
gain layer. The time-dependent radius and velocity of the shock can then be 
obtained as solutions of an initial value problem, which reflects the fact that the 
destiny of the shock depends on the initial conditions and is controlled by the 
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cumulative effects of neutrino energy deposition and mass accumulation in the 
gain layer Q. 

Such an analysis also demonstrates the existence of a threshold value for 
the neutrino luminosity from the neutron star, which is needed to drive shock 
expansion. This threshold luminosity depends on the rate, M, of mass infall to 
the shock, on the neutron star mass and radius, and to some degree also on the 
shock stagnation radius. Taking into account only the main dependence on M, 
it can roughly be written as |n| 

£„,crit(M) « Lo - L x (r^j > ( 13 ) 

with Lo w 5 x 10 52 ergs -1 and L\ w 3 x 10 52 ergs -1 for the conditions of Fig. |[ 
The neutrino heating in the gain layer is not the only important factor that 
determines the shock propagation. Energy loss by neutrino emission in the cool- 
ing layer has a considerable influence, because it regulates the settling of the 
matter that is accreted by the nascent neutron star, and therefore the advection 
of gas through the heating layer. If cooling is inefficient, gas piles up on the 
neutron star and pushes the shock farther out. If cooling is very efficient, the gas 
contracts quickly and more gas is dragged downward through the gain radius, 
extracting mass and energy from the gain layer and thus weakening the support 
for the shock. This also means that the infall velocity behind the shock increases 
and the timescale for the gas to stay in the gain layer is reduced. Therefore 
the efficiency of neutrino energy deposition drops. Such an effect is harmful for 
shock expansion. It can be diminished by higher v e and v e luminosities from the 
neutrinosphere, which lead to an enhancement of neutrino absorption relative to 
neutrino emission. On the other hand, muon and tau neutrino and antineutrino 
production in the accretion layer of the neutron star has a desastrous conse- 
quence for the shock, because it is a sink of energy that leaves the star without 
any significant positive effect above the neutrinosphere, where only v e and v e 
can be absorbed by free nucleons. 

These different, competing processes combined explain the slope of the criti- 
cal line in the left plot of Fig. g and in Eq. ( |l3| ) . Shock expansion and acceleration 
are easier for high mass infall rates, M, into the shock and for high v e and v e 
luminosities from the nascent neutron star. These luminosities need to be larger 
when M is small. It must be pointed out here, however, that this dependence is 
a consequence of the fact that the temperature in the cooling layer is considered 
as a parameter of the discussion. It is assumed to be equal to the neutrinospheric 
temperature and thus to be mainly determined by the interaction with the neu- 
trino flux from the core of the neutron star, but not by the mass infall and the 
dynamics in the accretion layer. 

Emission of muon and taun neutrinos and antineutrinos from the cooling 
region is not included in the results displayed in Fig. |[ It would move the critical 
line to appreciably higher values of the v e plus v e luminosity of the collapsed 
core, which is given along the ordinate. 
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Fig. 3. Explosion energies E>o{i) for spherically symmetric models ("ID", dashed 
lines) and two-dimensional ("2D") supernova models (solid lines) [^,^3[. The latter 
take into account convective overturn between the supernova shock and the neutrino- 
heating region. The curves display the evolution as a function of time after shock 
formation for different values of the v £ luminosity (labeled in units of 10 52 erg/s). The 
latter was used as a free parameter at the surface of the nascent neutron star and was 
roughly equal to the v e luminosity. Below the smallest given luminosities, the considered 
15 Mq star does not explode in ID and acquires too low an expansion energy in 2D to 
unbind the whole stellar mantle and envelope. 



Neutrino heating is stronger close to the gain radius than right behind the 
shock. Using an isentropic profile in the gain layer, the evaluation, however, 
implies very efficient energy transport, e.g., by convective motions in the gain 
layer. This enhances the postshock pressure and reduces the loss of energy from 
the gain layer, which is associated with the inward advection of neutrino-heated 
gas. 

Solutions of the toy model for varied parameters show that the energy in 
the gain layer and therefore the explosion energy of the supernova is limited to 
some 10 51 erg. The reason for this is the following. Neutrino energy deposition 
proceeds by v e and v e absorption on nucleons. The heated gas expands away 
from the region of strongest heating as soon as the nucleons have absorbed an 
energy roughly of the order of the gravitational potential energy, with only a 
small time lag because of the inertia of the shock, which is confined by the ram 
pressure of the collapsing stellar material. This does not allow the net energy of 
the heated gas to become very large. Typically it is of the order of ~ 5 MeV per 
nuclcon. With a mass in the gain layer of several 0.01 Mq up to ~ 0.1 Mq, the 
total energy does therefore not exceed a few 10 51 erg. 
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Figure q clearly shows this saturation of the explosion energy, which occurs 
when the gain layer and the shock are expanding. In this sense, neutrino-driven 
explosions are "self-regulated": Further energy deposition is quenched as the 
baryons move out of the region of high neutrino fluxes. 

The heating rate increases with the neutrino luminosity and the deposited 
energy is higher for a larger mass in the gain layer (Eq. ^). However, the ex- 
pansion timescale during which the gas is exposed to high neutrino fluxes, drops 
when the heating is stronger. Therefore the explosion energy is extremely sen- 
sitive to the neutrino luminosity only around the threshold value for getting an 
explosion. 

Neutrino-driven explosions are likely to be "delayed" (up to a few 100 ms 
after core bounce) rather than "late" (after a few seconds) . The density between 
the gain radius and the shock decreases with time because the proto-neutron 
star contracts and the gas infall to the shock drops rapidly as time goes on. 
Therefore the mass AM in the heated region shrinks and the shock must recede 
to a small radius around or even below 100 km, not favoring a later explosion. 

Fulfilling the "explosion criterion" of Fig. || once is no guarantee for a suc- 
cessful explosion. The push from the heating region has to be maintained until 
the material that carries the bulk of the energy is moving ballistically, and a 
fair fraction of this energy has been converted from internal to kinetic energy. 
Otherwise, if energy losses by neutrino emission or PdV on the contracting 
proto-neutron star start to dominate the energy input by neutrino absorption, 
the pressure-supported expansion can break down again and re-collapse can oc- 
cur. However, once shock expansion sets in, the conditions for further neutrino 
heating improve rapidly, and the optical depth of the growing gain layer to neu- 
trinos increases. Provided the neutrino luminosity does not drop, an explosion 
becomes unavoidable. This requirement favors a high core neutrino luminosity 
over accretion luminosity to power neutrino-driven explosions. 



3 Status of Spherical Simulations 

Recently, a major shortcoming of previous supernova models has been removed, 
at least in spherical models. Instead of treating multi-frequency neutrino trans- 
port by a flux- limited diffusion approximation []64|j65| , ^9| ,p| , the Boltzmann equa- 
tion can now be solved in connection with hydrodynamical simulations, either by 
direct discretisation |Q or by a variable Eddington factor technique |54j| , even 
in the general relativistic case fl^j . For the first time, the numerical deficiencies 
of the models are therefore smaller than the uncertainties of the input physics. 

The more accurate treatment of the transport, in particular in the semi- 
transparent neutrino-decoupling region around and outside of the neutrinosphere 
up to the shock, favors higher energy transfer to the stellar gas in the cooling 



layer and in the gain layer [|43 69 1. Nevertheless, the results of these simulations 
are disappointing. In spherical symmetry, the models do not explode, neither in 
the Newtonian (Fig. 0), nor in the general relativistic case. 
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Fig. 4. Trajectories of mass shells in the core of a collapsing 15 Mq star from a Newto- 
nian simulation with Boltzmann neutrino transport [p3|]54| l . The shells are equidistantly 
spaced in steps of 0.02 Mq. The boundaries of the iron core, silicon shell and neon- 
magnesium shell are indicated by bold lines. The fat, solid curve rising up at 0.21 
seconds after the start of the simulation marks the position of the supernova shock, 
the dashed line denotes the gain radius. In this spherically symmetric simulation, muon 
and tau neutrinos and antineutrinos were neglected, which favors an explosion. Never- 
theless, the shock recedes after having expanded to more than 300 km. 

These current models, however, neglect convective effects inside the nascent 
neutron star as well as in the neutrino-heated region. Since convection has been 
recognized to be very important, such simulations do not treat the full super- 
nova problem and do not really allow for conclusions about the viability of 
the neutrino-driven mechanism. Multi-dimensional simulations with Boltzmann 
neutrino transport are called for. 

4 Hydrodynamic Instabilities during the Explosion 

During the explosion of a supernova, hydrodynamic instabilities and convective 
processes can occur on different scales in space and time. Convective motions 
inside the nascent neutron star can speed up the energy transport and raise the 
neutrino luminosities during a period of seconds (Section 5). In the neutrino- 
heated region, convective overturn during the first second of the explosion car- 
ries hot matter towards the shock front and brings cool gas into the region of 
strongest neutrino heating near the gain radius. This has important influence on 
the start of the explosion and the nucleosynthesis of radioactive elements. When 
the shock propagates through the mantle and envelope of the disrupted star, 
Rayleigh- Taylor instabilities destroy the onion-shell structure of the progenitor 
and mix radioactive material with high velocities from near the neutron star 
into the helium and even hydrogen shells of the star. In this section, the early 
postshock convection and its interaction with the hydrodynamic instabilities at 
the composition interfaces of the progenitor star will be discussed. 
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Fig. 5. Inhomogeneous distribution of neutrino-heated, hot gas which rises in 
mushroom-like bubbles, and cooler gas that is accreted through the supernova shock 
(bumpy discontinuity at about 3000 km) and falls in long, narrows streams towards the 
newly formed neutron star at the center. The figure shows a snapshot of the entropy 
at 300 ms after core bounce and shock formation for a two-dimensional simulation of 
a 15 Mq star. The star was exploded by the neutrino-heating mechanism by chosing 
a suitable value of the neutrino luminosity from the nascent neutron star, which was 
replaced by an inner boundary condition ^(]^| . The white line encompasses the region 
where radioactive nickel has been formed by nuclear burning in the shock-heated Si 
layer. 



4.1 Convective Overturn in the Neutrino-Heated Region 

Convective instabilities in the layers adjacent to the nascent neutron star are a 
natural consequence of the negative entropy gradient built up by the weakening 
of the prompt shock prior to its stagnation and by neutrino heating Ml. This was 
verified by two- and three-dimensional simulations 0,|^,|8|j2^,|8|,|8|,||j3|Q . 
Figure |^ shows the entropy distribution between proto-ncutron star and super- 
nova shock about 300 ms after core bounce for one such calculation |4C|^l[ |. 
Although there is general agreement about the existence and the growth of hy- 
drodynamic instabilities in the layer between the shock at R s and the radius 
of maximum neutrino heating (which is just outside the gain radius, R g ), the 
strength of the convective overturn and its importance for the success of the 
neutrino-heating mechanism are still a matter of debate. 

Two-dimensional simulations with a spectrally averaged, flux-limited diffu- 
sion treatment of neutrino transport |28| , |20[] , or with the neutrino luminosity be- 
ing given as a free parameter at the inner boundary, which replaces the neutron 
star at the center [ ^2||33]| , found successful explosions in cases where spherically 
symmetric models fail (Fig. |^). According to these simulations, the convective 
overturn in the neutrino-heated region has the following effects on the shock 
propagation. Heated matter from the region close to the gain radius rises out- 
ward and at the same time is exchanged with cool gas flowing down from the 
shock. Since the production reactions of neutrinos {e^ capture on nucleons and 
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thermal processes) are very temperature sensitive, the expansion and cooling 
of rising plasma reduces the energy loss by the reemission of neutrinos. More- 
over, the net energy deposition by neutrinos is enhanced as more cool material 
is exposed to the large neutrino fluxes near the gain radius (the radial dilution 
of the fluxes goes roughly as l/r 2 ). Since hot matter moves towards the shock, 
the pressure behind the shock increases, an effect which pushes the shock farther 
out. This leads to a growth of the gain region and therefore also of the net energy 
transfer from neutrinos to the stellar gas, favoring an explosion. 

The consequences of postshock convection are clearly visible from the re- 
sults plotted in Fig. ||, where the explosion energy E >Q as a function of time is 
shown for spherically symmetric and two-dimensional calculations of the same 
post-collapse model, but with different assumed neutrino luminosities from the 
proto-neutron star |||,|33| ■ ^>o is defined to include the sum of internal, kinetic, 
and gravitational energy for all zones where this sum is positive (the gravitational 
binding energies of stellar mantle and envelope and additional energy release 
from nuclear burning are not taken into account). For one-dimensional simula- 
tions with v e luminosities (and very similar v e luminosities) below 1.9- 10 52 erg/s 
explosions could not be obtained when the proto-neutron star was assumed 
static, and the threshold value of the luminosity was 2.2 • 10 52 erg/s when the 
neutron star was contracting. The supporting effects of convective overturn be- 
tween the gain radius and the shock lead to explosions even below the critical 
value in spherical symmetry, and to a faster development of the explosion. 

Simulations with a better description of the neutrino transport by a multi- 
energy-group treatment of the neutrino diffusion pE[ , confirm the existence of 
such convective processes in the region of neutrino heating, but the associated 
effects are not strong enough to revive the stalled prompt supernova shock, 
although the outward motion of the shock is enhanced. 

Fully self-consistent, multi-dimensional calculations, however, have not yet 
been done with a state-of-the-art Boltzmann neutrino transport, which has re- 
cently become applicable for spherically symmetric models (see Section 3). The 
current multi-dimensional simulations therefore demonstrate only the presence 
and potential importance of convection, but final conclusions on the viability of 
the neutrino-heating mechanism in the presence of postshock convection are not 
possible at the moment. A quantitatively meaningful description of the shock 
revival phase, however, requires an accurate description of the transport as well 
as a multi-dimensional approach. 

4.2 Is Neutrino-Driven Convection Crucial for an Explosion? 

The role of convective overturn for the development of an explosion becomes 
clearer by considering the three timescales of neutrino heating, r^t, advection 
of accreted matter through the gain layer into the cooling region and down to 
the neutron star (compare Fig. Q), r a d, and the timescale for the growth of 
convective instabilities, r cv . The evolution of the shock — accretion or explosion 
- is determined by the relative size of these three timescales. Straightforward 
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Fig. 6. Order scheme for the post-collapse dynamics in dependence of L v {e^}, which 
determines the strength of the neutrino heating outside of the neutrinosphere. The 
destiny of the star — accretion or explosion — can be understood by the relative size 
of the timescales of neutrino heating, Tht, matter advection through the gain region 
onto the nascent neutron star, r a d v , and growth of convective instabilities, r cv . 



considerations show that they are of the same order and the destiny of the star 
is therefore a result of a tight competition between the different processes. 

The heating timescale is estimated from the initial entropy per nucleon, sj, 
the critical entropy s c (Eq. pQ)), and the heating rate per nucleon (Eq. (0)) as 



s c - Si Sc - sj J% 7 (T/2MeV) / 

™ ~ Qt/(k B T) ~ DmS ' 5k B /N (L,/4-10^erg/s)(e2 15 ) ' [ > 

for L v being the total luminosity of v e plus v e . With a postshock velocity of 
til = uq/ ' [3 ~ (7 — l)y/GM/Rs/ (7 + 1) the advection timescale is 

Rst — R e /„ Re\ ^s,200 



r a d « « 55 ms • 1 - ^ , (15) 

mi \ R S J 



where the gain radius can be determined as 



-1/4 
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from the requirement that the heating rate, Eq. (|4|), is equal to the cooling rate 
per nucleoli, Q~ k, 2.3 T 6 MeV/(N ■ s), when -R s ,200 is the shock radius in units 
of 200 km, Y n + Y p w 1 is assumed, and use is made of the power-law behavior 
of the temperature according to T(r) « T s (R s /r), with T s being the postshock 
temperature in MeV. The growth timescale of convective instabilities in the 
neutrino-heated region depends on the gradients of entropy and lepton number 
through the growth rate of Ledoux convection, ctl : 



dp\ ds ( dp \ dY c 



dsjy^ pAr \dY e J sP dr 



-i/2 



- p \ 1/2 R 3/2 

20ms- (#-1) _M_ _ (17) 



^g 

The numerical value was obtained with the gravitational acceleration g = GMj R^, 
(dp/ds)p ~ —p/s, and ds/dr ~ — is/(i? s — i? g ). The term proportional to the 
gradient of Y" e was assumed to be negligible. t cv of Eq. ([l?]) is sensitive to the 
detailed conditions between gain radius (close to which s develops a maximum) 
and the shock. The neutrino heating timescale is shorter for larger values of 
the neutrino luminosity L v and mean squared neutrino energy (e 2 ,). All three 
timescales, Tht, T a d and r cv , decrease roughly in the same way with smaller gain 
radius or shock position. 

In order to be a crucial help for the explosion, convective overturn in the 
neutrino-heated region must develop on a sufficiently short timescale. This hap- 
pens only in a rather narrow window of L„(ef,) where r cv < r a d ~ Tht (Fig. ^). 
For smaller neutrino luminosities the heating is too weak to create a sufficiently 
large entropy maximum, and rapid convective motions cannot develop before 
the accreted gas is advected through the gain radius (r a d < "?cv < "Hit)- In this 
case neither with nor without convective processes energetic explosions can oc- 
cur (Fig. [5j). For larger neutrino luminosities the neutrino heating is so strong, 
and the heating timescale correspondingly short (Tht < r cv ~ 7"ad)i that expan- 
sion of the postshock layers has set in before the convective activity reaches a 
significant level. In this case convective overturn is an unavoidable side-effect 
of the neutrino heating behind the shock, but is not necessary for starting the 
explosion. 

The parametric studies performed by Janka and Miiller fl32|]33|l support this 
discussion, which helps one understanding the seemingly discrepant results ob- 
tained by different groups. 



4.3 Nucleosynthesis and Mixing Instabilities 

Besides increasing the efficiency of neutrino energy deposition, convection in the 
postshock layer has an important influence also on the nucleosynthesis and dis- 
tribution of radioactive elements. In particular, nickel is not produced by silicon 
burning in a spherical shell, but is concentrated in dense clumps and pockets 
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Fig. 7. Snapshot from a two-dimensional simulation of the explosion of a blue su- 
pergiant star with 15 solar masses at a time 1170 seconds after the stellar core has 
collapsed to a neutron star In the left half of the figure the density is shown in a 
region with a radius of about 2.2 million kilometers, in the right half three color im- 
ages of the mass densities of radioactive nickel (red and pink), silicon (green, light blue, 
whitish) and oxygen (deep blue) are superposed. One can see that the ejecta of the 
explosion are inhomogeneous and anisotropic, and the original onion-shell structure of 
the exploding star was shredded. Nickel is concentrated in dense, fast-moving clumps 
along the extended filaments seen in the left plot. 




Fig. 8. Left: Distribution of Ni vs. radial velocity at t — 50 s and t = 300 s after 
core bounce. Right: Initial composition of the star exterior to the iron core (top) and 
composition 300 s after core bounce (bottom). C, O, Si, and the newly synthesized Ni 
have been mixed beyond the inner half of the helium core, and He has been carried 
inward. 



between rising bubbles of neutrino-heated matter in the expanding postshock 
layer (Fig. |). 

The further evolution of the shock until it breaks out of the stellar sur- 
face hours later, was recently followed by using adaptive mesh refinement tech- 
niques These allow for a dynamic adjustment of the computational grid 
such that small structures can be treated with high resolution, while the whole 
computation covers a huge volume. 

A few seconds after its formation, the shock has passed the silicon and 
oxygen layers and propagates through the helium shell of the star. The initial 
anisotropies have been compressed into a narrow shell, from which new instabil- 
ities start to grow. Rayleigh- Taylor mushrooms penetrate into the helium layer 



18 H.-Th. Janka, K. Kifonidis, and M. Rampp 



and carry 0, Si and Ni farther out, while He sinks in. Within minutes, long, dense 
filaments reach far into the helium shell, associated with them fast-moving knots 
that contain dominant contributions of different heavy elements from the deeper 
layers (Fig. [?]). Nickel, silicon and oxygen move through helium with velocities 
up to several 1000 km/s (Fig. ||). 

These simulations show that the hydrodynamic instabilities which occur in 
the first second of the explosion, do not only play a role during the phases of 
shock rejuvenation and nickel formation. They act also as seed perturbations 
for the instabilities at the composition interfaces of the progenitor, which finally 
destroy the onion-shell structure of the pre-collapse star. 

The results are in good agreement with observations of mixing and anisotropics 
in many Type Ib,c supernovae. In case of Supernova 1987A, a Type II explosion 
of a massive star which has retained its hydrogen envelope, the observed high 
nickel velocities in the hydrogen envelope cannot be explained by the models. 
The nickel clumps are strongly decelerated at the He/H interface, where they 
enter a dense helium "wall" which builds up after the passage of the shock. The 
dissipation of the kinetic energy of the clumps does not allow nickel to penetrate 
into the hydrogen layer with high velocities. 



5 Neutron Star Formation 

Convective energy transport inside the newly formed neutron star can increase 
the neutrino luminosities considerably | ft7j . This can be crucial for energizing 
the stalled supernova shock []64|65f| . 

Convection in the neutron star can be driven by gradients of the entropy 



and/or proton (electron lepton number) fraction in the nuclear medium 22 . The 
type of instability which grows most rapidly, e.g., doubly diffusive neutron- finger 
convection or Ledoux convection [jl7j or quasi-Ledoux convection (3^^§] , 

may be a matter of the properties of the nuclear equation of state, which deter- 
mines the magnitudes and signs of the thermodynamic derivatives [[l3| . It is also 
sensitive to the gradients that develop, and thus may depend on the details of 
the treatment of neutrino transport in the dense interior of the star. 

Convection below the neutrinosphcre seems to be disfavored during the very 
early post-bounce evolution by the currently most elaborate supernova mod- 
els 1 14 1^,0 , but can develop deeper inside the nascent neutron star on a longer 
timescale (kJ 100 ms after bounce) and can encompass the whole star within 
seconds [ p7|j3^ , |36| . 

Negative lepton number and entropy gradients have been seen in several one- 
dimensional (spherically symmetric) simulations of the neutrino-cooling phase 
of nascent neutron stars Jl6|,[l7 37,5^] (see also Fig. ^) and have suggested the 
existence of regions which are potentially unstable against Ledoux convection. 
Recent calculations |50| with improved neutrino opacities of the nuclear medium, 
which were described consistently with the employed equation of state, confirm 
principal aspects of previous simulations, in particular the existence of Ledoux- 
unstable layers in the neutron star. 
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Fig. 9. Profiles of the lepton fraction Yj ep = ni ep /nt (left) and of the entropy per 
nucleon, s, (right) as functions of enclosed (baryonic) mass for different times in a 
one-dimensional simulation of the neutrino cooling of a ~ 1.1 Mq proto- neutron star. 
Negative gradients of lepton number and entropy suggest potentially convectively un- 
stable regions. Time is (roughly) measured from core bounce. 



5.1 Convection inside the Nascent Neutron Star 

Two-dimensional, hydrodynamical simulations were performed for the neutrino- 
cooling phase of a ~ 1.1 Mq proto- neutron star that formed in the core collapse of 
a 15 Mq star p8pr| . The models followed the evolution for a period of more than 

1.2 seconds. They demonstrate the development of convection and its importance 
for the cooling and deleptonization of the neutron star. 

The simulations were carried out with the hydrodynamics code Prometheus. 
A general relativistic ID gravitational potential with Newtonian corrections for 
asphcricities was used, <P = + ($2D ~ ^id)> anc ^ a flux- limited (equilibrium) 
neutrino diffusion scheme was applied for each angular bin separately ("l^D"). 

The simulations showed that convectively unstable surface-near regions (i.e., 
around the neutrinosphcre and below an initial density of about 10 12 g/cm 3 ) 
exist only for a short period of a few ten milliseconds after bounce, in agreement 
with the findings by other groups fl4| , |l5l , |44f . Due to a flat entropy profile and 
a negative lepton number gradient, convection, however, also starts in a layer 
deeper inside the star, between an enclosed mass of 0.7 M© and 0.9 Mq, at 
densities above several 10 12 g/cm 3 . From there the convective region digs into 
the star and reaches the center after about one second (Figs. 0, and [Hf). 
Convective velocities as high as 5-10 8 cm/s were found (about 10-20% of the local 
sound speed), corresponding to kinetic energies of up to 1-2 • 10 50 erg (Fig. |TT| ) . 
Because of these high velocities and rather flat entropy and composition profiles 
in the star (Fig. |l(J), the overshooting region is large (see Fig. |T3| ). The same is 
true for undershooting during the first ~ 100 ms after bounce. Sound waves and 
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Fig. 10. Same as Fig. Q but for a two-dimensional, hydrodynamical simulation which 
allowed to follow the development of convection. The plots show angularly averaged 
quantities in the ~ 1.1 Mq proto-neutron star. In regions with convective activity the 
gradients of Y\ ep and s are flattened. The convective layer encompasses an increasingly 
larger part of the star. 



Fig. 11. Absolute values of the convective velocity in the proto-neutron star for two 
instants (about 0.5 s (left) and 1 s (right) after core bounce) as obtained in a two- 
dimensional, hydrodynamical simulation. The arrows indicate the direction of the ve- 
locity field. Note that the neutron star has contracted from a radius of about 60 km 
initially to little more than 20 km. The growth of the convective region can be seen. 
Typical velocities of the convective motions are several 10 s cm/s. 



perturbations are generated in the layers above and interior to the convection 
zone. 

The coherence lengths of convective structures are of the order of 20-40 
degrees (in 2D!) (see Fig. [ll]) and coherence times are around 10 ms, which 
corresponds to only one or two overturns. The convective pattern is therefore 
very time-dependent and nonstationary. Convective motions lead to considerable 
variations of the composition. The lepton fraction (and thus the abundance of 
protons) shows relative fluctuations of several 10%. The entropy differences in 
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Fig. 12. Convective "luminosity" (solid line) and neutrino luminosities (dashed: L„ e + 
L^ e , dash-dotted: L„ M +-^p h + L Vt + Lq t , dotted: total) as functions of enclosed baryonic 
mass for the two-dimensional proto-neutron star simulation about 500 ms after core 
bounce. 

rising and sinking convective bubbles are much smaller, only a few per cent, 
while temperature and density fluctuations are typically less than one per cent. 

The energy transport in the neutron star is dominated by neutrino diffusion 
near the center, whereas convective transport plays the major role in a thick 
intermediate layer where the convective activity is strongest. Radiative transport 
takes over again when the neutrino mean free path becomes large near the surface 
of the star (Fig. [l2]). But even in the convective layer the convective energy flux 
is only a few times larger than the diffusive flux. This means that neutrino 
diffusion is not negligibly small in the convective region. This fact has important 
consequences for the driving mechanism of the convection. 

5.2 Driving Force of Convection 

The convective activity in the neutron star cannot be explained by, and consid- 
ered as ideal Ledoux convection. Applying the Ledoux criterion for local insta- 
bility, 

with o"l from Eq. (|l7|) and Y e replaced by the total lepton fraction Y\ cp in the 
neutrino-opaque interior of the neutron star (for reasons of simplicity, Vs was 
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Fig. 13. Left: Convectively unstable region (corresponding to negative values of the 
displayed quantity ljq L = — (<?/p)Cql with Cql from Eq. (|l9|)) about 500 ms af- 
ter bounce according to the Quasi-Ledoux criterion which includes non-adiabatic and 
lepton-transport effects by neutrino diffusion. Right: Layer of Quasi-Ledoux convective 
instability (blue) as function of time for a two-dimensional simulation. The angle- 
averaged criterion Cql(r) = mine (Cql (r, 6)) > with Cql (r, 0) from Eq. ( p^| ) is 
plotted. The dotted area is outside of the computed star, green denotes stable layers 
where over- and undershooting causes lateral velocities with angularly averaged abso- 
lute values of (\vg\) > 10 7 cms _1 , and white are convectively "quiet" regions of the 
star. 

replaced by ds/dr and VFi ep by dY\ cp /dr), one finds that the convecting re- 
gion should actually be stable, despite of slightly negative entropy and lepton 
number gradients. In fact, below a critical value of the lepton fraction (e.g., 
Yj epc = 0.148 for p = 10 13 g/cm 3 and T = 10.7 MeV) the thermodynamical 
derivative [dp/ dY\ cp ) S: p changes sign and becomes positive because of nuclear 
and Coulomb forces in the high-density equation of state [OJ . Therefore nega- 
tive lepton number gradients should stabilize against convection in this regime. 
However, an idealized assumption of Ledoux convection is not fulfilled in the 
situations considered here: Because of neutrino diffusion, energy exchange and, 
in particular, lepton number exchange between convective elements and their 
surroundings are not negligible. Taking the neutrino transport effects on Yj ep 
into account in a modified "Quasi-Ledoux criterion" |3(| , 

(19) 

one determines instability exactly where the two-dimensional simulation reveals 
convective activity. In Eq. ( |l9|) the quantities (Y\ ep ) and (s) mean averages over 
the polar angles 6, and local gradients have to be distinguished from gradients 
of angle-averaged quantities which describe the stellar background. The term 
Aep(diiep/dr) with the empirically determined value f3\ ep ~ 1 accounts for the 
change of the lepton concentration along the path of a rising fluid element due 
to neutrino diffusion. Figure |l3] shows that about half a second after core bounce 
strong, driving forces for convection occur in a narrow ring between 9 and 10 km, 
where a steep negative gradient of the lepton fraction exists (see Fig. [h]). Farther 
out, convective instability is detected only in finger-like structures of rising, high- 
Yi ep gas. 
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Fig. 14. Absolute value of the gas velocity in a convecting, rotating proto-neutron 
star about 750ms after bounce (left). Convection is suppressed near the rotation axis 
(vertical) and develops strongly only near the equatorial plane where a flat distribution 
of the specific angular momentum j z (right) has formed. 



5.3 Accretion and Rotation 

In other two-dimensional models, post-bounce mass accretion and rotation of 
the forming neutron star were included. Accretion causes stronger convection 
with larger velocities in a more extended region. This can be explained by the 
steepening of lepton number and entropy gradients and the increase of the grav- 
itational potential energy when additional matter is added onto the neutron 
star. 

Rotation has very interesting consequences, e.g., leads to a suppression of 
convective motions near the rotation axis because of a stabilizing stratification of 
the specific angular momentum (see Fig. |l4|), an effect which can be understood 
by applying the (first) Solberg-H0iland criterion for instabilities in rotating, self- 
gravitating bodies 



dx 



•st) I 
1 djf , a 



5s /y- lcp ,P S \dY lcp J sP lcp 



< . (20) 



Here, j z is the specific angular momentum of a fluid element, which is conserved 
for axially symmetric configurations, x is the distance from the rotation axis, and 
in case of rotational equilibrium a is the sum of gravitational and centrifugal 
accelerations, a = VP/ p. Changes of the lepton number in rising or sinking 
convective elements due to neutrino diffusion were neglected in Eq. (|2^). Ledoux 
(or Quasi-Ledoux) convection can only develop where the first term is not too 
positive. In Fig. |lj fully developed convective motion is therefore constrained to 
a zone of nearly constant j z close to the equatorial plane. At higher latitude the 
convective velocities are much smaller, and narrow, elongated convective cells 
aligned with cylindrical regions of j z = const parallel to the rotation axis are 
visible. 

The rotation pattern displayed in Fig. |lj is highly differential with a rotation 
period of 7.3 ms at x — 22 km and of 1.6 ms at x = 0.6 km. It has self-consistently 
developed under the influence of neutrino transport and convection when the 
neutron star had contracted from an initial radius of about 60 km (with a surface 
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Fig. 15. Left: Luminosities L v (i) and mean energies of u e and v e for a 1.1 Mq 

proto-neutron star without ("ID"; dotted) and with convection ("2D"; solid). Right: 
Angular variations of the neutrino flux at different times for the 2D simulation. 

rotation period of 55 ms at the equator and a rotation period of ~ 5 ms near the 
center) to a final radius of approximately 22 km. Due to the differential nature of 
the rotation, the ratio of rotational kinetic energy to the gravitational potential 
energy of the star is only 0.78% in the beginning and a few per cent at the end 
after about 1 s of evolution. 

5.4 Consequences of Proto-Neutron Star Convection 

Convection inside the proto-neutron star can raise the neutrino luminosities 
within a few hundred ms after core bounce (Fig. |l5|). In the considered collapsed 
core of a 15 M@ star, L Vc and Lp e increase by up to 50% and the mean neutrino 
energies by about 15% at times later than 200-300 ms post bounce. This favors 
neutrino-driven explosions on timescales of a few hundred milliseconds after 
shock formation. Also, the deleptonization of the nascent neutron star is strongly 
accelerated, raising the v e luminosities relative to the v e luminosities during this 
time. This helps to increase the electron fraction Y e in the neutrino-heated ejecta 
and might solve the overproduction problem of N = 50 nuclei during the early 
epochs of the explosion fl38|] . In case of rotation, the effects of convection on 
the neutrino emission depend on the direction. Since strong convection occurs 
only close to the equatorial plane, the neutrino fluxes are convectively enhanced 
there, while they are essentially unchanged near the poles. 

Anisotropic mass motions due to convection in the neutron star lead to grav- 
itational wave emission and anisotropic radiation of neutrinos. The angular vari- 
ations of the neutrino flux found in the 2D simulations are of the order of 5-10% 
(Fig. [if]). With the typical size of the convective cells and the short coherence 
times of the convective structures, the global anisotropy of the neutrino emission 
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Fig. 16. Kick velocity of the neutron star as a function of time, caused by the 
anisotropic emission of neutrinos due to convection. The two-dimensional simulation 
was done with a polar grid from to n. 

from the cooling proto-neutron star is very small. This implies a kick velocity of 
the nascent neutron star due to anisotropic neutrino emission of only ~ 10 km/s 
in a 2D simulation (Fig. |l^). Because the convective elements are likely to be- 
come even smaller in 3D, kick velocities of 300 km/s or even more, as observed for 
many pulsars, can definitely not be explained by convectively perturbed neutrino 
emission. 

5.5 Neutrino Opacities in Nuclear Matter and Neutron Star 
Convection 

Another important issue of interest are the neutrino opacities in the dense and 
hot nuclear medium of the nascent neutron star. In current supernova models, the 
description of neutrino-nucleon interactions is incomplete because the standard 
approximations assume isolated and infinitely massive nucleons [ |S2"| . Therefore 
effects like the fermion phase space blocking of the nucleons, the reduction of the 
effective nucleon mass by momentum-dependent nuclear interactions in the dense 
plasma, and nucleon thermal motions and recoil are either neglected completely 
or approximated in a more or less reliable manner ]9|,^6|. These effects have 
been recognized to be important ]57| , p^ , p r l|j55[ | for calculations of the neutrino 
luminosities and spectra, but still await careful inclusion in supernova codes. For 
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Fig. 17. Left: Thermal averages of the neutrino mean free paths for absorption 
(dotted line), v<, scattering (solid line), and muon and tau neutrino scattering (dashed 
line), respectively, according to the standard description of the neutrino opacities. M 
is the baryonic mass enclosed by the radial coordinate R (bold solid line) of the newly 
formed neutron star about 200 ms after core bounce. Also shown are the factors 1 and 
/y 1 (dash-dotted lines) which give a measure of the increase of the neutrino mean free 
paths caused by a suppression of the axial- vector and vector current contributions to the 
neutrino opacities due to in-medium effects (Yamada, personal communication). 
Right: Total neutrino luminosities, L v , and integrated energy loss, E v , as functions of 
time for spherically symmetric models without convection (Bl and B1YO) and two- 
dimensional models with convection (B2 and B2YO). Models Bl and B2 were computed 
with standard neutrino opacities whereas in B1YO and B2YO in-medium suppression 
of the neutrino opacities was included [feq] (Yamada, personal communication) . 



this purpose a consistent description of nuclear equation of state and neutrino- 
matter interactions is desirable. 

Many-body (spatial) correlations due to strong interactions []56|,^9 19[55| , |68f 
and multiple-scattering effects by spin-dependent forces between nucleons (tem- 
poral spin-density correlations) |52 26 1 are of particular interest, because they 
lead to a reduction of the neutrino opacities in the newly formed neutron star 
and are associated with additional modes of energy transfer between neutrinos 
and the nuclear medium. 

A reduction of the neutrino opacities implies larger neutrino mean free paths 
and thus increases the neutrino luminosities. (Fig. and Refs. j3§|^||]). The 
neutrino diffusion is accelerated most strongly in the very dense core of the 
nascent neutron star. Convection in the intermediate region between core and 
outer layers turns out not to be suppressed, but is still the fastest mode of 
energy transport. Therefore reduced neutrino opacities as well as convective 
energy transport are important, but the combined effects do not appreciably 
change the convectively enhanced neutrino emission (Fig. [l^, right) |3q] . 
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6 Summary 

Supernova explosions of massive stars are an important phenomenon for apply- 
ing nuclear and particle physics, in particular neutrino physics. The processes 
going on in the extremely dense and hot core of the exploding star are accessible 
to direct measurements only through neutrinos or gravitational waves. Empir- 
ical information about the events that cause the explosion and accompany the 
formation of a neutron star, however, can also be deduced from observable char- 
acteristics of supcrnovae, for example their explosion energy or the amount and 
distribution of radioactive nuclei, and from the properties of neutron stars. 

Theoretical models need to establish the link between the core physics and 
these observables. In the past years it has been recognized that hydrodynam- 
ical instabilities and mixing processes on large scales play an important role 
within the core as well as in the outer layers of the exploding star. Convection 
can change the cooling of the nascent neutron star, supports the revival of the 
stalled shock by neutrino heating, and destroys the onion-shell structure of the 
progenitor star. Multi-dimensional calculations are therefore necessary to un- 
derstand why and how supcrnovae explode, and to make predictions for their 
observable consequences. 

Spherically symmetric simulations, Newtonian and general relativistic, with 
the most advanced treatment of neutrino transport by solving the Boltzmann 
equation, do not produce explosions. This emphasizes the importance of convec- 
tion, but may also point to physics still missing in the models. One such weakness 
of current simulations is an overly simplified description of neutrino interactions 
with nuclcons in the nuclear medium of the neutron star. A kinematically cor- 
rect treatment of these reactions, taking into account nucleon thermal motions, 
recoil and fermi blocking, needs only a technical step, but a better understand- 
ing of the effects of nucleon correlations and their consistent treatment with the 
equation of state requires theoretical progress. 

The neutrino-heating mechanism, although the favored explanation for the 
explosion, is still controversial, both because of the status of modeling and be- 
cause of observations which seem hard to explain. Although significant progress 
has been made, multi-dimensional simulations with an accurate and reliable han- 
dling of neutrino transport and an up-to-date treatment of the input physics are 
still missing, and definite conclusions can therefore not be drawn at the moment. 
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